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Abstract 



[The low-lying collective states in Sn isotopes are studied by a five-dimensional collective Hamiltonian with parameters determined 

from the triaxial relativistic mean-field calculations using the PC-PK1 energy density functional. The systematics for both the 

O-iexcitation energies of 2t states and B(E2; Ot — > 21") values are reproduced rather well, in particular, the enhanced E2 transitions 
flj ■ i ii 

in the neutron-deficient Sn isotopes with N < 66. We show that the gradual degeneracy of neutron levels lgj/2 and 2d$/2 around 

the Fermi surface leads to the increase of level density and consequently the enhanced paring correlations from N = 66 to 58. 

.It provokes a large quadrupole shape fluctuation around the spherical shape, and leads to an enhanced collectivity in the isotopes 

^ aroundJV = 58. 

' ' Keywords: covariant energy density functional, collective Hamiltonian, low-lying states, electromagnetic transition, Sn isotopes 
+2 IpACS: 21.60.Jz, 21.60.Ev, 21.10.Re, 21.10.Tg 
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The evolution of shell structure and nuclear collectivity along 
either isotopic or isotonic chains, inferred from the system- 
atics of the nuclear deformation energy surface HI], the nu- 
cleon separation energy 

a a, 

the excitation energy of low- 
lying state and the electric quadrupole transition strength 

an, 

becomes of great interest in nuclear physics. In recent years, 
a number of experiments have measured the B(E2; 0| — > 2|) 
.transition strengths in neutron-deficient Sn isotopes with N < 
66 fi0 HSCl! El El ■ The deduced B(E2) value increases 



when going from N - 66to N - 64, and then remains roughly 
constant within the experimental uncertainties when decreasing 
the neutron number down to N = 56. This picture deviates 
from the parabolic trend, with a peak at mid-shell, predicted by 
the single y'-shell seniority model [13]. Furthermore, the B(E2) 
values deduced from the measurement of the lifetimes of the 
first excited 2 + states in I12 . 114 .H6sn at the UNILAC accelerator 
of the Gesellschaft fur Schwerionenforschung (GSI) are overall 
smaller than the previously reported values and present a shal- 
low minimum at N = 66 Il4ll . These findings have generated a 
renewal interest in the study of low-lying states in Sn isotopes. 

The theoretical description of the enhanced B(E2; 0[ — > 2|) 
transition strength in the neutron-deficient Sn isotopes is not 
straightforward. Similarly as the seniority scheme, the senior- 
ity truncated large-scale shell model calculations predicted a 
nearly parabolic behavior of B(E2) values with a maximum at 
N = 66 • By using new interactions and including core po- 
larization terms up to both the third order and Shut core ex- 
citations in the shell model calculations, the predicted B(E2) 
values exhibit a shallow maximum around N = 68 - 70, but 



are still far away from the experimental data Holl . Shell model 
calculations in the sdgh major shell using CD-Bonn and Ni- 
jmegenl two-body effective nucleon-nucleon (NN) interactions 
were also carried out for the mid-heavy even-even Sn isotopes 
m Refs. An overall good agreement with the observed 

energy spectra of Sn isotopes is found. However, the E2 transi- 
tions are not presented. Recently, a quasiparticle-phonon model 
(QPM) was adopted to investigate the evolution of the 2| state 
in Sn isotopes 11711 . The calculation reproduced the trend of the 
energies and, partly, the observed deviations of the E2 strengths 
from the parabolic behavior. It was shown that such an asym- 
metric trend was the consequence of several factors: single- 
particle energies, polarization of the N - Z - 5Q core, interplay 
between pairing plus quadrupole, and quadrupole pairing inter- 
actions II 1711 . all of which, however, are quite model-dependent. 

In the framework of nuclear energy density functional theory 
or self-consistent Hartree-Fock-Bogoliubov (HFB) approach, 
due to the robust Z = 50 shell gap, the ground states of all the Sn 
isotopes are dominated by the spherical configuration. For the 
2\ states, which are expected to be a one-phonon state caused 
by surface vibration around the spherical shape, the theoret- 
ical studies are obviously beyond the mean-field approxima- 
tion. For this purpose, one can either introduce random-phase 
approximation (RPA) or generator coordinate method (GCM), 
where the polarization of N = Z = 50 core in the shell model 
can be taken into account by particle-hole excitation in RPA or 
configuration mixing of different shapes in GCM. The imple- 
mentations of the quasiparticle RPA on top of the nonrelativis- 
tic and relativistic mean-field calculations were carried out to 
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investigate the evolution of low-lying states along the Sn iso- 
topic chain 0, [Hill 0] . In Refs. Ql^, the systematics 
of B(E2; 0| — > 2|) values can be roughly reproduced. How- 
ever, these studies were focused on the substantial rise of the 
B(E2) values at neutron magic numbers 50 and 82, and did not 
provide any interpretation for the enhanced collectivity in the 
neutron-deficient Sn isotopes. 

In recent years, the GCM has been introduced into the self- 
consistent mean-field approaches by combing with projection 
techniques in the modern energy density functional calcula- 
tions lilHQESIUliilzl]. The dynamic correla- 



tion effects related to the symmetry restoration and quadrupole 
fluctuation (along both B and y directions) around the mean- 
field minimum could be taken into account. These methods 
restricted to axial case have been applied to study the low- 
lying collective excitation spectra of the neutron-deficient Pb 
isotopes ifjol 31]. However, the application of these meth- 
ods with triaxial degree-of-freedom for systematic study is still 
much time-consuming. Up to now, such kind of study is mostly 
restricted to l igh t nuclei 11321 13311 and some specific medium 
heavy nuclei 11341. 13511 . As the Gaussian overlap approximation 
of GCM, collective Hamiltonian with parameters determined 
by self-consistent mean-field calculations is much simple in nu- 
merical calculations, and has achieved great success in descrip- 
tion of nuclear low-lying states for deformed and transitional 
nuclei Ezl USE? 



41,42,43,44,45 



In the past decades, the covariant density functional theory 
(CDFT) has achieved great success in the description of nu- 
clear structure over almost the whole nuclide chart, from light 
systems to superheavy nuclei, and from the valley of B stability 
to the particle drip lines |47, 48, 4^]. Recently, the implemen- 
tation of the five-dimensional collective Hamiltonian (5DCH) 
based on the CDFT has been developed and given successful 
microscopic description for the nuclear sharje transitions along 
both the isotopic and isotonic chains f5l [38, 4(| 41, 42]. In 
this work, we are going to apply this approach to study the 
low-lying collective states, for the first time, in the semi-magic 
Sn isotopes. Our aim is to provide a microscopic mechanism 
for the enhanced E2 transitions in the neutron-deficient Sn iso- 
topes. 

In the calculation, we use a recent parameterized relativis- 



tic functional PC-PK1 [50] for the particle-hole channel, and 
a separable force 15 ll 15211 for the particle-particle channel. 



The solution of the equation of motion for the nucleons is ac- 
complished by an expansion of the Dirac spinors in a set of 
three-dimensional harmonic oscillator basis functions in Carte- 
sian coordinates with 12 major shells. More details about the 
calculations can be found in Ref. 1271 l38l l53f . The intrin- 



sic triaxially deformed states are obtained as solutions of the 
self-consistent RMF+BCS equations constrained on the mass 
quadrupole moments related to the Bohr parameters (B, y) vary- 
ing B e [0.0, 0.6] and y e [0°, 60°] with step size AB = 0.05 and 
Ay = 10°. 

The nuclear collective excitations associated with three- 
dimensional rotation and quadrupole fluctuations are described 
with a five-dimensional collective Hamiltonian (5DCH). The 
dynamics of the 5DCH is governed by the seven functions of the 



intrinsic deformations B and y: the collective potential V co ii, the 
three mass parameters: Bpp, Bp r , B yy , and the three moments of 
inertia 1^. These functions are determined using cranking ap- 
proximation formula based on the intrinsic triaxially deformed 
states. The diagonalization of the Hamiltonian yields the ex- 
citation energies and collective wave functions that are used to 
calculate observables 13811. 
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Figure 1: (Color online) Calculated £(2|) (upper panel) and B(£"2;0| 
(lower panel) in Sn isotopic chain by the shell model with 100 Sn core (dashed 
line), RQRPA (dotted line), and 5DCH based on PC-PK1 density functional 
(solid line with open circles), in comparison with the data (solid symbols). In 
the lower panel, the measured B(E2\0^ —> 2|) values from different experi- 
ments are indicated with different symbols. 

The upper panel of Fig. [TJ displays the excitation energies 
of 2| states for Sn isotopes from our 5DCH calculations, in 
comparison with those by the shell model with 100 Sn core 
and RQRPA calculations, as well as the data. The trend of 
the excitation energies of 2| states along the isotopic chain 
is reproduced quite well by the three models. It is notable 
that our 5DCH calculations give quite similar results as those 
by the RQRPA calculations, although the realization of these 
two models for excitation states are quite different. Compared 
with the shell model calculation, the discrepancy between our 
5DCH result and the data is relatively larger for 120 ~ 130 Sn, from 
~ 200 keV to ~ 1 MeV. However, it is worthwhile to note that 
the 5DCH is fully self-consistent and without any free param- 
eter. In addition, we find that, beside 2| state, the excitation 
energies of other low-lying states for 108-1 16 Sn are described 
as well as that by the recent shell model calculations lfl5l [l6Tl . 
Taking 116 Sn as an example one finds that the rms deviation for 
seven low-lying excitation states (the first two + , three 2 + , and 
two 4 + ) is: cr = 0.279 MeV for the 5DCH, which is in compari- 
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son with the corresponding value 0.429 MeV obtained by using 
Nijmegenl effective interaction in Ref. [15], and 0. 137 MeV in 
the calculation of Ref. lfl6Tl . 

In the lower panel of Fig. Q] we present the B(E2; 0| — > 2|) 
values from 5DCH calculations, in comparison with the results 
by shell model and RQRPA, as well as the recently measured 
data. It is shown that both the 5DCH and RQRPA calculations 
reproduce rather well the systematics of B(E2; 0[ — > 2|) val- 
ues in the whole isotopic chain, in particular the increase of E2 
strength from I16 Sn down to 108 Sn, which cannot be described 
by the shell model. Quantitatively, the RQRPA slightly under- 
estimates those E2 strengths around I16 Sn. While the 5DCH 
overestimates the E2 strength ~ 0.1 e 2 b 2 when approaching to 
the neutron shell closure, where the low-lying states are mainly 
caused by the excitation of few quasiparticles. Concerning our 
interest and the power of our model, in the following we will 
mainly focus on the isotopes around middle shell, where the ex- 
citation is dominated by the collective vibration and enhanced 
E2 transition is observed. 
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of neutron number. A rapid drop of the curvature in Sn isotopes 
with A < 116 is found, which is responsible for the increased 
B(E2) value shown in Fig. Q] To understand this change in the 
topograph of deformation energy curves, we plot the average 
neutron pairing gaps, weighted by uv coefficient as Ref. |54[], 
at spherical states in the lower panel of Fig. [3] For compari- 
son, the odd-even mass differences of experimental data by the 
five-point formula are also plotted. It is shown that the rapid 
drop of the curvature when decreasing neutron number from 
N = 66 is closely correlated to the neutron pairing gaps. Simi- 
lar evolution trend is shown in the odd-even mass differences. It 
indicates that the strong pairings soften the energy curve around 
spherical shape, which leads to the enhanced collectivity in Sn 
isotopes around N = 58. 
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Figure 2: (Color online) Self-consistent RMF+BCS binding energy curves for 
the even-even 108-1 16 Sn isotopes as functions of the axial deformation parame- 
ter fi. Energies are normalized with respect to the binding energy of the corre- 
sponding spherical state. 



To understand the enhanced E2 transitions around N = 
108, we plot the deformation energy curves of the even-even 
108-1 i6g n i so t p es as functions of the axial deformation param- 
eter B in Fig. [2] Due to the robust Z = 50 proton shell, all the 
Sn isotopes are spherical or near spherical. Around the neutron 
mid-shell N = 66, an excited prolate minimum with B ~ 0.3 
shows up. As the decrease of the neutron number from N - 66 
down to N = 58, the energy curve around spherical shape be- 
comes softer gradually. In particular, for 108 Sn, a weakly prolate 
deformed global minimum with B ~ 0.1 is observed. In other 
words, the dynamic effect of quadrupole deformation fluctua- 
tions is enlarging the nuclear collectivity gradually when de- 
creasing neutron number from N - 66 down to N - 58. 

The upper panel of Fig.[3]displays the curvatures of the defor- 
mation energy curves around the spherical shape as a function 



Figure 3: (Color online) The predicted curvature of the PEC (panel a) and av- 
erage neutron pairing gaps (panel b) at spherical point in Sn isotopic chain 
calculated by PC-PK1 density functional. For comparison, the odd-even mass 
differences (denoted by solid squares) extracted by the five-point formula are 
also plotted in panel b. 

Figure |4] displays the neutron single-particle energies at 
spherical shape along Sn isotopic chain. The balls indicate 
the position of Fermi surface. It is seen that the levels lg-j/2 
and 2c/5/2 move closer as the neutron number decreases from 
N = 82 to N = 50 shell. Meanwhile, the Fermi surface 
comes toward to these two levels for the isotopes with A < 116. 
This leads to the increasing level density around Fermi surface 
and consequently enhanced paring correlations. Moreover, it is 
found that the Fermi surface crosses the 3si/2 orbit in between 
N - 6A and N = 70. The considerable occupation probability 
of 3ii/2 orbit hinders the collectivity of these nuclei. Mean- 
while, the obstructive effect is further enhanced by the low den- 
sity of single-particle states around the Fermi surface, which 
reduces the pair correlations (c.f. in Fig. [3]). It provides a qual- 
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Figure 4: (Color online) Evolution of neutron single-particle levels at spherical 
point for Sn isotopic chain. The balls indicate the position of Fermi surface. 



itative interpretation of the dip of E2 transition at 1 16 Sn, which 
is consistent with that given in Ref . IU4I1 . 

In summary, the low-lying collective states in Sn isotopes 
from N = 50 to N = 82 have been studied by solving a 
five-dimensional collective Hamiltonian with parameters deter- 
mined from the triaxial relativistic mean-field calculations. The 
calculations reproduce the systematics of both excitation en- 
ergies for 2\ states and B(E2; 0[ — > 2|) rather well. In par- 
ticular, the enhancement of E2 transitions from N = 66 to 
N = 58 has been found to be the consequence of increasing 
quadrupole shape fluctuation around spherical shape, induced 
by the strengthening neutron pairing correlations. This con- 
clusion is consistent with that by the shell models, but in the 
language of mean-field approach. Furthermore, it is found the 
gradual degeneracy of levels lgj/2 and 2^5/2 around the Fermi 
surface that leads to the increase of level density and conse- 
quently the enhanced paring correlations. Moreover, it has 
shown that the Fermi surface crosses the 3si/2 orbit in between 
N = 64 and N = 70, where the low level density and consid- 
erable occupation probability of 3s\/2 orbit hinder the nuclear 
collectivity and give rise to valley-like of the E2 transition as a 
function of neutron number. 
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